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We found 140 neutrino-induced muons in 854.24 live days in the MINOS far detector, which has 
an acceptance for neutrino-induced muons of 6.91 x 10 6 cm 2 sr. We looked for evidence of neutrino 
disappearance in this data set by computing the ratio of the number of low momentum muons to the 
sum of the number of high momentum and unknown momentum muons for both data and Monte 
Carlo expectation in the absence of neutrino oscillations. The ratio of data and Monte Carlo ratios, 
1Z, is 

71 = 0.65t°;l2(stat) ± 0.09(syst), 

a result that is consistent with an oscillation signal. A fit to the data for the oscillation parameters 
sin 2 2#23 and Am| excludes the null oscillation hypothesis at the 94% confidence level. We separated 
the muons into /i~ and /x + in both the data and Monte Carlo events and found the ratio of the total 
number of to /i + in both samples. The ratio of those ratios, TZcpt, is a test of CPT conservation. 
The result 

Kcpt = 0.72«; 2 *(stat)t°;^(syst), 
is consistent with CPT conservation. 

PACS numbers: 



I. INTRODUCTION 

Measurements of atmospheric neutrinos by Super- 
Kamiokande experiment have shown that there is a 
deficit of when compared to expectations p], 0, 0, 0, 0, 
B 0, B 0] ■ The hypothesis that best describes this deficit 
is the oscillation of v^v^) — > v T (p T ) [HI HH, with the 
oscillation probability given by 

P v ^ Vt = sin 2 26» 23 sin 2 (1.27 Am 2 23 L/E), (1) 

where 023 is the mixing angle, Am 2 3 = |m 2 — m 2 | is 
the mass squared difference in eV 2 between the neu- 
trino mass states, L is the distance in km traveled by 
the neutrino, or its baseline, and E is the energy of 
the in GeV [T3|. The Super-Kamiokande data is 
best fit by the oscillation hypothesis with parameters 
(sin 2 26»23,Am| 3 ) = (1.0, 2.4x 10" 3 eV 2 ) @, Q ; the ranges 
for these parameters given by the 90% confidence con- 
tours of the zenith angle oscillation fit are sin 2 2023 > 0.92 



*Now at Cleveland Clinic, Cleveland, Ohio 44195, USA. 
tNow at GE Healthcare, Florence South Carolina 29501, USA. 
'Now at APC — College de France, 11 Place Marcelin Berthelot, 
F-75231 Paris Cedex 05, France. 

§Now at Physics Department, Royal Holloway, University of Lon- 
don, Egham, Surrey, TW20 0EX, United Kingdom. 
^Now at Holy Cross College, Notre Dame, Indiana 46556, USA. 
"Now at Lincoln Laboratory, Massachusetts Institute of Technol- 
ogy, Lexington, Massachusetts 02420, USA. 
"Deceased. 

M Now at Nuclear Nonproliferation Division, Threat Reduction Di- 
rectorate, Los Alamos National Laboratory, Los Alamos, New Mex- 
ico 87545, USA. 



and 1.5 < Am^^K 3.4 x 10" 3 eV 2 . The MACRO [H03 
and Soudan 2 [15|, results are consistent with those 
obtained by Super-Kamiokande. For the MINOS anal- 
ysis of atmospheric neutrinos with an interaction vertex 
in the detector, the parameter ranges are sin 2 2023 > 0.2 
and 7 x 10" 5 < Am| 3 < 5 x 10~ 2 eV 2 Below 
we extend the MINOS atmospheric analysis to neutrino- 
induced muons observed in the MINOS far detector. 

The oscillation hypothesis for the atmospheric neu- 
trino deficit has received strong support from the first 
results of the MINOS long baseline experiment. MI- 
NOS, which sends produced at Fermilab to a de- 
tector 735 km away in northern Minnesota (HI, finds 
(sin 2 26»23, Amy = (1.00, 2.74 xl0~ 3 eV 2 ) [H, when fit- 
ting their events; the ranges of these parameters given 
by the 90% confidence contours are sin 2 2023 > 0.72 and 
2.2 < Am 2 3 < 3.8xlO" 3 eV 2 . The K2K long baseline ex- 
periment, which sent produced at KEK to the Super- 
Kamiokande detector located 250 km away, also found 
consistent results - (1.0, 2.8 x 10" 3 eV 2 ) [H El[ with 90% 
confidence ranges for the parameters of sin 2 20 2 3 > 0.55 
and 1.9 < Am| < 3.6 x 10~ 3 eV 2 . 

Since the magnetized MINOS far detector distin- 
guishes from MINOS data can be used as a probe 
of CPT conservation in the neutrino sector. CPT con- 
servation requires that the probability of an atmospheric 
of a given energy remaining a after traveling from 
its origin to its point of interaction be equal to the prob- 
ability that an atmospheric 77^ of the same e nerg y re- 
mains a 7/ M after traveling the same distance [ll|. We 
consequently test CPT conservation by comparing v^- 
induced \iT and F p -induced /i + with respect to expecta- 
tions, as a measure of whether the atmospheric and 
v ^ are disappearing at the same rate. Effects induced by 
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charged current interactions of the v e components of the 
neutrino eigenstates with the earth could masquerade as 
apparent CPT violation. However, distortions from mat- 
ter effects are estimated to be small for most of the #13 
range compatible with the CHOOZ limit ■ The aver- 
age induced v^jv^ event rate asymmetry is estimated to 
be less than 10% [23J. Therefore, at our current level of 
statistical accuracy we can ignore matter effects and de- 
scribe v^iv^) — » v T (p T ) oscillations using the two flavor 
approximation in Eq. (TTJ). 

After briefly discussing the detector in SjTTl we describe 
our event sample in ^III Al and the Monte Carlo event 
generation in Bl In ^IIICl we describe the cuts 
used in selecting the sample of muons to be analyzed. 
We present the oscillation analysis in §IVI and the probe 
of CPT conservation using charge separated neutrino- 
induced fx~~ and /i + in i fVl 



4000 




II. 



THE MINOS FAR DETECTOR 



2 2.5 

RMS (ns) 

FIG. 1: Distribution of the r.m.s. deviations of the measured 
hit times from the calculated hit times for calibration tracks. 
The resolution of the system is the mean of a Gaussian fit to 
the distribution, 2.31 ± 0.03 ns. 



The MINOS far detector is a steel-scintillator sampling 
calorimeter located at a depth of 2070 meters-water- 
equivalent (m.w.e.) in the Soudan mine in northern Min- 
nesota [l8j. The detector is made up of 486 vertical 
octagonal planes of 2.54 cm thick steel laminates, inter- 
leaved with 484 planes of 1 cm thick extruded polystyrene 
scintillator strips and a 2.5 cm wide air gap. Each scintil- 
lator plane has 192 strips of width 4.1 cm. The length of 
each strip depends on its position in the plane and varies 
between 3.4 and 8.0 m. The scintillator strips in alter- 
nating detector planes are oriented at ±45° to the verti- 
cal. The modular detector consists of two supermodules 
(SM) separated by a gap of 1.1 m. The detector coordi- 
nate system has the y axis pointing up, the z axis along 
the detector long axis, pointing away from Fermilab and 
the x axis forms a right-handed coordinate system. 

Scintillation light from charged particles traversing the 
MINOS plastic scintillator is collected with wavelength 
shifting (WLS) plastic fibers embedded in the scintillator 
strips. The WLS fibers are coupled to clear optical fibers 
at both ends of a strip and are read out using 16-pixel 
multi- anode phototmultiplier tubes (PMTs). The signals 
from eight strips, separated by approximately 1 m within 
the same plane, are optically summed, or multiplexed, 
and read out by a single PMT pixel. The multiplexing 
pattern is different for the two sides of the detector, en- 
abling the resulting eightfold ambiguity to be resolved for 
single particles. For all other types of events, ambiguities 
are resolved effectively using additional information from 
timing and event topology. 

To measure the momentum of muons traversing the 
detector, the steel has been magnetized into a toroidal 
field configuration using a current-carrying coil that runs 
through the central axis of each SM. A finite element 
analysis calculation shows that each SM is magnetized to 
an average value of ~ 1.3 T by the 15 kA current loop. 
The field is saturated near the coil hole at a strength of 



approximately 1.8 T, falling to about 1 T near the edges. 
Further details about the MINOS detector can be found 




To distinguish upward-going neutrino-induced muons 
from the downward cosmic ray background requires ex- 
cellent timing resolution. The MINOS timing system has 
a least count of 1.56 ns. However, the timing resolution 
is dominated by the 8 ns fluorescence decay time in the 
WLS fiber 24]. The overall resolution is determined by 
the convolution of the least count with the fluctuations 
in the arrival times of the photons at the PMT. We cal- 
ibrated the timing system by measuring the time offsets 
between each channel. For this determination [25|, [26| . 
we selected a sample of downward through-going cosmic 
ray muons with well-reconstructed tracks. For each hit 
along a track we measured the travel times from the track 
entrance point in the detector to the hit location and cor- 
rected those values for both the rise time of the signal 
and the propagation of light along the fibers. We then 
compared these times with the time expected for a rela- 
tivistic muon traversing the same distance and computed 
the offset between the measured and expected times for 
each hit, At. The timing calibration constants are ob- 
tained using an iterative method to minimize the mean 
At for each channel. After calibration, linear timing fits 
are again applied to the times and positions of the hits 
on each muon track. The r.m.s. deviation between the 
measured and fitted times is calculated for each track. 
Figure [T] shows the distribution of these r.m.s. deviations. 
In this figure, the distribution has been fit to a Gaussian 
with a mean of 2.31 ± 0.03 ns, a value that represents the 
overall timing resolution of the detector. 
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III. 



DATA ANALYSIS 



In the analyses presented here we use neutrino-induced 
muons, which are defined as events that come from below 
or slightly above the horizon. These events are essen- 
tially uncontaminated by the background of downward- 
going atmospheric muons. We analyze two types of 
events: those that pass completely through the detec- 
tor ("through-going muons") and those that enter and 
stop in the detector ("stopping muons"). Earlier MINOS 
results for events with an interaction vertex in the detec- 
tor ("contained events") appear in [T7| . We distinguish 
between neutrino-induced muons with well-determined 
momenta and neutrino- induced muons whose charge and 
momentum are undetermined (c.f., Dp . 



A. Data Sample 

The data presented here represent 854.24 live days and 
were taken between August 1, 2003 and April 30, 2006. 
The geometric acceptance of the detector for neutrino- 
induced muons is 6.91 x 10 6 cm 2 sr. We found a total of 
140 neutrino-induced muons in this data set. 



B. Simulated Muons 




u.\ Super-K 

Oscillation Parameters 

H*, Super-K 

Oscillation Parameters 



(GeV/c) 



FIG. 2: Input momentum distribution for neutrino-induced 
yT and fi + calculated with the Nuance simulation package. 
The distributions in the MINOS Far detector without oscil- 
lations and with the Super-Kamiokande determined values of 
sin 2 2#23 and Am| 3 are shown. The oscillation affects mostly 
muons with low momenta. 



to see whether it survived energy loss processes between 
the surface and the detector [31| . We placed the muons 
that survived on an imaginary box surrounding the detec- 
tor and then propagated them though the detector 
with the MINOS GEANT-based detector simulation. 



1. Atmospheric Neutrino-induced Muons 

We generated a large sample of simulated neutrino- 
induced muons in the MINOS far detector using the 
Bartol 96 [13] neutrino flux tables along with the de- 
fault Nuance neutrino interaction model [28| and the 
Grv94 [2^] parton distribution functions. We first sim- 
ulated neutrino-induced muons passing through the sur- 
face of a box surrounding the MINOS far detector and 
then propagated these muons through the detector with 
the MINOS GEANT-based detector Monte Carlo simula- 
tion. A total of 6.5 x 10 5 neutrino-induced muon events, 
the equivalent of 2500 years of live time, were generated 
in this way. The momentum distributions simulated for 
the neutrino-induced [i~ and /i + in MINOS are shown in 
Fig. [2] for neutrinos without oscillations and for neutrinos 
with an oscillation signal using the Super-Kamiokande 
parameters [8]. This figure shows that the oscillation 
signal affects mostly muons with low momenta. 



2. Cosmic Ray Muons 

We simulated each cosmic ray muon event by first 
choosing its arrival direction from a distribution uni- 
form in solid angle and then associating this direction 
with the overburden found in the Soudan 2 slant depth 
map (3(ij . The surface energy of the muon was selected 
from the known distribution [3l| . Once the surface en- 
ergy and overburden were chosen, we tested the muon 



C. Event Selection 

All data and simulated muons were analyzed using the 
standard MINOS reconstruction algorithms and a uni- 
form set of event selection cuts. The selection cuts sum- 
marized below are described in further detail in 12511. 



1. Muon Selection Cuts 

We first selected muon events using criteria developed 
for the study of cosmic ray muons with the MINOS far 
detector 



The first two cuts require that there was a 
track successfully reconstructed in the event ( "No Track" 
cut), and that there is only a single track found by the 
track-fitting algorithm ("Multiples" cut). The next set 
of cuts exclude random collections of hits that could be 
mistaken for muon tracks. These cuts require that: the 
track must cross at least 20 planes in the detector ("20 
Plane" cut); the track must have a pathlength of at least 
2 m ("2.0 m Length" cut); the earliest recorded hit of 
the track be no more than 15 cm from the front or back 
surface of the detector or 30 cm from any of the remaining 
surfaces ("Fiducial" cut); and the track fit must have a 
X 2 per number of degrees of freedom, x 2 fit /ndf < 1.5. 

The next selection cuts were designed to remove events 
with poor timing information. To ensure the presence of 
sufficient timing information in the events, we first ex- 
cluded tracks if fewer than half of the hits come from 
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strips with signals on both ends ("Double-ended Strip 
Cut"). We determined the direction of the track by plot- 
ting the time difference AT (ns) of each hit along the 
track as a function of its distance AS (m) from the first 
hit. If the y positions of the hits increase along the length 
of the track, AS is positive; for y decreasing along the 
track, AS is negative. The slope of the linear fit to the 
AT/ AS distribution is 1/(3 = c/v. The top panel of 
Fig. [3] shows the AT/ AS distribution with the linear fit 
superposed for a typical cosmic muon; the bottom panel 
displays the (a;, y) positions of the hits in the detector for 
this muon. 

In Fig. [4] we show the distributions of x^^/ndf values 
for the linear fits to 1/(3 for both data and cosmic ray 
muon Monte Carlo events passing the above cuts. The 
distribution does not peak at 1 as might be expected 
for a x 1 /ndf statistic. However, that is unimportant 
to the analysis since we are using the x\/p/ n< V value 
only to define a cut. We defined a cut that requires 
xl/p/ndf < 3.0 for an event ("xup/ndf < 3.0" cut). 

The cut was selected at the value of Xi/p/ndf where the 
number of events falls to ~1% of the peak. This cut 
maximizes the number of selected events while excluding 
those events with the worst Xi/p/ndf values. Although 
the data and Monte Carlo simulation deviate at low val- 
ues of x\/p/ n df the agreement is excellent near the peak 
and the cut value. Therefore, the Monte Carlo sample 
can be used to effectively study the systematic uncer- 
tainty introduced into the analysis by this cut. 

In the last step of the muon selection, we separated 
the downward-going from upward-going events. Upward- 
going events have a positive slope for the straight line fit 
to the AT / AS distribution. The final muon cut is a 
check on the up/down separation. The entrance point 
of a track, and therefore its incoming direction, is deter- 
mined by the slope of the times of the hits in the detector 
as a function of their z position. For a few short tracks 
this information is sufficiently ambiguous that the en- 
trance point can be confused with the exit point of the 
track causing the reconstruction to interchange the two. 
As a result, downward-going muons are incorrectly re- 
constructed as upward-going. We remove these failures 
by checking the times of the hits in the reconstructed 
track as a function of their y position. The slope of the 
line fit to the hits as a function of their y positions must 
agree with the reconstructed incoming direction of the 
track ("Directionality" cut). 

Table U shows the effect of the cuts on the data and 
the two Monte Carlo distributions. We normalize the 
cosmic ray muon Monte Carlo distributions to the data 
value at the 20 Plane cut because the Monte Carlo sim- 
ulation does not account for multiple muon events which 
make up ~ 5% of downward-going cosmic ray muons [34| . 
After cuts the selected events consist of muon tracks 
whose directionality is well determined. The distribution 
of 1/(3 values for the selected muon events is shown in 
Fig. [H Since the sample includes both downward-going 
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FIG. 3: A typical cosmic ray muon in the MINOS far detector. 
The top panel shows the timing information for the hits along 
the track with a straight line fit superposed. The legend gives 
AT /AS = 1//3 and the xl/fj/ndf of the fit for this muon. The 
bottom panel shows the {x, y) hit positions of this track. The 
resolution of the positions is the 4 cm width of the strips. 



and upward-going muons and these muons are relativis- 
tic, there are two peaks in the 1 / (3 distribution. The peak 
at 1/(3 = — 1 are downward-going muons and the peak at 
1/(3 = 1 are upward-going muons. Integrating the events 
in the two peaks shows that the fraction of (upward-going 
muons)/(downward-going muons) ~ 10 -5 . 



2. Upward-Going Neutrino-Induced Muons 

As in previous experiments [3, fl(| , we identify upward- 
going muons that enter the detector as neutrino-induced 
muons. Events in the range 0.7 < 1/(3 < 1.3 are included 
in our neutrino-induced muon data set. The cut was 
determined using downward-going cosmic ray muons. It 
was set as the 1/(3 range that includes 99% of the events 
and which is centered on the peak. In the left panel of 
Fig. [S] we show a well-characterized upward-going muon 
event. The right panel shows an event excluded by the 
1/(3 cut. 

Table |TT] shows that our data sample includes 130 
upward-going, neutrino-induced muons in the range 
0.7 < 1/(3 < 1.3. In this table we also show the back- 



FIG. 4: Distribution of Xi/p/ndf values from the 1/(3 deter- 
minations for data and cosmic ray muon Monte Carlo events. 
In this figure, the Monte Carlo distribution has been normal- 
ized to the same area as the data distribution. Events with 
Xi/p/ndf < 3 are used in the analysis. 



TABLE I: Effect of cuts on data and simulated samples. The 
fraction of the total events remaining is shown. N gives the 
total number of events in each sample. 



Cut 


Data 


/u MC 


i/-induced MC 




N = 


N = 


N = 




3.81 x 10 7 


1.20 x 10 6 


6.47 x 10 s 


No Cuts 


1.000 




1.000 


No Track 


0.800 




0.862 


Multiples 


0.743 




0.848 


20 Plane 


0.561 


0.561 


0.606 


2.0 m Length 


0.557 


0.557 


0.578 


Fiducial 


0.534 


0.538 


0.559 


Xm/ndf < 1.5 


0.429 


0.447 


0.497 


Double-ended Strips 


0.429 


0.447 


0.497 


Xt/p/ndf < 3.0 


0.428 


0.447 


0.497 


Directionality 


0.428 


0.447 


0.478 
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FIG. 5: Distribution of 1//3 for upward-going neutrino- 
induced muons, with a peak at 1, and downward-going cosmic 
ray muons, with a peak at -1. The vertical lines at 0.7 and 
1.3 bracket the events included in the upward-going muon 
sample. 



FIG. 6: Distributions of AT/ AS as in Fig. El The left panels 
shows a typical upward-going muon and the right panel shows 
a muon excluded by the 1//3 cut from the upward- going sam- 
ple. 



ground from contained vertex interactions where the neu- 
trino interacts close to the detector edge and cannot 
be separated from muons entering the detector. This 
background was estimated by calculating the rate of 
Monte Carlo atmospheric neutrino events that interact 
inside the detector volume [171 ]. The background in- 
cludes events that pass the cuts in Table Q] and satisfy the 
upward-going selection cut. The Monte Carlo events were 
assumed to oscillate with the best-fit oscillation param- 
eters from Super-Kamiokande Q. The 20 Plane cut has 
been shown to be effective at removing the background 
due to upward-going pions produced by downward-going 
cosmic ray muons interacting in the rock surrounding the 
detector. 



3. Horizontal Muons 

In this analysis, we also include muons coming from 
along and slightly above the detector's horizon. The flat 
overburden of the Soudan site [30j makes this search 
feasible. The slant depth of rock between the detec- 
tor and the surface for incoming directions above the 
horizon increases approximately as sec#, where 9 is the 
zenith angle. Only the highest energy cosmic ray muons 
(Efj, > 100 TeV) have sufficient energy to penetrate the 
large column of rock present for muons coming from di- 
rections near the horizon. Since the intensity of cosmic 
ray muons falls as E~ 2 ' 7 , very few cosmic ray muons sur- 
vive to reach the detector, implying that muons from 
near the detector's horizon are neutrino-induced. These 
events are important to the analysis because they sample 
neutrino- induced muons with lower values of L/E than 
the upward- going muons. Eq. ([T]) shows that the average 
oscillation probability for these horizontal muons is lower 
than for upward-going neutrino-induced muons and they 
add information important in determining the neutrino 
flux normalization, a parameter used in the oscillation 
analysis in ^IV CI 
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TABLE II: Observed upward-going and horizontal muons. 
The background accounts for both contained vertex neutrino 
interactions and downward-going cosmic ray muons. 



Muon Type | Events 


Background 




Contained 


Cosmic /i 


Upward-going 


130 


4.2 


0.0 


Horizontal 


10 


0.1 


0.3 




FIG. 7: The cos# distribution for all incoming muons. The 
data are fit to an exponential between 0.1 < cos# < 0.2 to es- 
timate the background of cosmic ray muons in the horizontal 
neutrino-induced muon signal region. The fit is shown in the 
figure. As indicated, we select muons coming from directions 
with cos 6 < 0.05. 



We require that the measured velocity of the 
downward-going muons be in the range —1.3 < 1//3 < 
—0.7 . To determine the maximum angle above the hori- 
zon from which we can distinguish horizontal muons from 
cosmic ray muons, we use the zenith angle distribution 
for all events shown in Fig. [71 The distribution is steeply 
falling for cos# < 0.25 but becomes approximately con- 
stant for cos# < 0.1. The rate of neutrino-induced muons 
is expected to be approximately independent of cos 9 [3l| . 
We select events with cos (9 < 0.05 to minimize the back- 
ground from cosmic ray muons. To estimate this back- 
ground we fit an exponential to the distribution of Fig. [7J 
for 0.1 < cos# < 0.2, a similar procedure to that used 
by Super-Kamiokande The exponential fit is shown 
in the figure; it has the form y — exp(a -I- bcosO), where 
a = —3.59 and b = 51.65. Fig. [7] shows that the back- 
ground at cos 9 — 0.05 is 0.3 events and is negligible for 
smaller values of cos 9. 

In Table |TT] we show that there arc 10 horizontal, 
neutrino-induced muons in the data sample. The back- 
ground shown includes both cosmic ray downward-going 
muons and the contained vertex events discussed previ- 
ously. 



D. Charge Sign and Momentum Determination 

Previous oscillation analyses based on neutrino- 
induced muons have typically divided the event sample 
into through-going and stopping muons. The experi- 
ments that presented these analyses could only measure 
the muon momentum directly for stopping muons. Since 
the MINOS far detector is magnetized, we can measure 
the muon momentum and charge sign for both types of 
neutrino-induced muons. 

To determine whether the charge and momentum of a 
neutrino-induced muon has been determined accurately, 
we have developed a criterion based on the measured 
curvature of the muon track. We first draw a straight 
line between the end points on the track and then use the 
deviations of the remaining track points from the line to 
form a x 2 statistic, X^ ine /ndf. Large values of 'xl inc /ndf 
indicate significant bending in the magnetic field which 
leads to a good determination of the charge sign and 
momentum. 

We use our Monte Carlo simulation to set the value of 
X\ine/ n< V that selects neutrino-induced muons with good 
momentum and charge sign determination. As a test of 
whether the X\i n e/ n( V values are well simulated we com- 
pare in Fig. [5] the distributions of x\nJ n< V f° r stopping 
muons in the cosmic ray data, in the cosmic ray muon 
Monte Carlo simulation, and in the neutrino-induced 
muon Monte Carlo simulation. The Monte Carlo distri- 
butions have been normalized to have the same number 
of events as in the data. Fig.[5]shows that the cosmic ray 
muon Monte Carlo simulation provides a high statistics 
match to the cosmic ray data. Since the stopping cosmic 
ray muons have energies < 20 GeV, the maximum energy 
of a normally incident stopping muon, these low energy 
cosmic events provide a sample of muons with a distribu- 
tion in energy similar to neutrino- induced muons. Also 
superposed on these distributions is that of the Monte 
Carlo neutrino-induced muons. This figure suggests that 
the neutrino-induced Monte Carlo simulation can be used 
to determine an appropriate value of x^ ne /ndf for select- 
ing events with good charge sign and momentum deter- 
mination. 

In Fig. [9] we show the fraction of Monte Carlo neutrino- 
induced muons with correct charge identification, or pu- 
rity, as a function of X 2 ine /nci/. The purity is approx- 
imately 97% at xjjne/Wf = 10 and rises to over 99% 
with increasing values of X}me/ n( V- We identify muons 
with X}me/ nc tf > 10 as having good charge-sign and mo- 
mentum determination; muons with X\incl n< V < 10 are 
assumed to have poor charge-sign and momentum deter- 
mination. 

We further divided the Monte Carlo muons with 
X\\nel n df > 10 into two samples based on their mo- 
mentum, a separation that roughly distinguishes muons 
whose parent neutrinos have a relatively large proba- 
bility of oscillation from those with a lower probability. 
The 'low momentum' (L) sample includes muons with fit 
momentum, Pfu, in the range 1 < put < 10 GeV/c. 
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FIG. 8: Comparison of the Xiine/ n df distributions for stop- 
ping rnuons in the cosmic ray data (points) , in the cosmic ray 
Monte Carlo simulation (boxes), and in the neutrino-induced 
Monte Carlo simulation (line). The vertical extent of the 
boxes representing the cosmic ray Monte Carlo simulation in- 
dicate the statistical uncertainty for the points. 
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FIG. 9: Purity of the charge sign determination for neutrino- 
induced muons as a function of y$\ uc /ndf. Events with 
Xiine/ n df > 10 have well-determined charge sign and momen- 
tum. 



The 'high momentum' (H) sample includes muons in 
the range 10 < Pfu < 100 GeV/c. No muons with 
xfi ne /ndf > 10 have pf it > 100 GeV/c. The muons with 
Xhne/ n( V < 10 are in the 'unknown momentum' (U) sam- 
ple. 

The quality of the momentum determination can 
be deduced from the Monte Carlo simulation us- 
ing the r.m.s. momentum resolution, (Ap/p) = 

y/((ptruth ~ Pfit) 2 /Ptruth)' wnere Ptruth is the known 
momentum of the muons. In Fig. [10] we show Ap/p, 
as a function of Xii n e/ nc tf f° r the low momentum and 
high momentum samples. The low momentum Monte 
Carlo muons have a momentum resolution of less than 
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FIG. 10: Momentum resolution for Monte Carlo neutrino- 
induced muons as a function of Xiinc/ n df- The line at 
Xiine/ n df = 10 shows the cut used in the analyses. 



50% at X\inel n df — 10 an d the resolution decreases 
to 10% with increasing values of x^ ino /nri/. For the 
high momentum muons the momentum resolution is less 
than 30% at Xi 2 ine /«c{f = 10 and it decreases to 10% as 
XimJndf increases. 

Fig. [TT] shows the distributions of parent neu- 
trino energies for the low momentum, high momentum, 
and unknown momentum muons as determined by the 
Monte Carlo simulation. The neutrinos producing 
low momentum muons have energies that peak near 
10 GeV. The high momentum muons are produced by 
neutrinos with energies that peak near 50 GeV. The 
unknown momentum muons are produced by even higher 
energy neutrinos, with a peak energy near 250 GeV. 
Comparing the peak energies of the various neutrino dis- 
tributions we see from Eq. ([T]) that the low momentum 
muons are those that are expected to show the strongest 
oscillation signal, a result also seen in Fig. [2] 

In Fig. [12] we show the distributions of x^ inc /nrif for 
the selected neutrino-induced muons from the data and 
unoscillated Monte Carlo simulation normalized by live 
time. The muons with Xii n e/ nc V < 10 arc from high 
energy neutrinos and are not expected to show a strong 
oscillation signal. Those muons with X 2 ine/ n( V > 10 are 
from lower energy neutrinos, a sample expected to show 
an oscillation signal. 

The results of applying the X\inol n ^f cu t to the sam- 
ple of 140 neutrino-induced muons are given in Table [TTT1 
In this table, the events with Xiine/ nc tf > 10 have been 
further separated by charge identification and momen- 
tum. The calculated background contribution is shown 
as well as the Monte Carlo expectation in the absence of 
oscillations. Figure [T31 shows the fit momentum distribu- 
tions for data and unoscillated Monte Carlo simulation 
using the combined low momentum and high momentum 
muon samples. The unknown momentum muons are not 
included in this figure. The first bin in Fig. [T3] shows 
a deficit of events in the data relative to Monte Carlo 
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FIG. 11: Distribution of energies for neutrinos producing 
neutrino-induced muons observed in the MINOS detector as 
determined by the Monte Carlo simulation. The neutrinos 
producing low momentum muons are shown by the solid line, 
those producing high momentum muons are shown by the 
dashed line and those producing unknown momentum muons 
are shown by open circles. Neutrinos producing muons de- 
tected by MINOS have energies > 2 GeV. 
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FIG. 12: Comparison of the x}\ ue /ndf distributions for 
neutrino-induced muon data and unoscillated Monte Carlo 
simulation. 



expectation without oscillations (c.f., Fig. 



IV. OSCILLATION ANALYSIS 



After first discussing the systematic uncertainties asso- 
ciated with the analysis procedure, we test the neutrino- 
induced muons in Table Mil for evidence of neutrino os- 
cillations. 



FIG. 13: Distribution of fit momenta for events in the com- 
bined low momentum and high momentum data samples. 
The Monte Carlo expectation for no oscillations is shown by 
the solid line. The unknown momentum muons are not in- 
cluded in this figure. 



TABLE III: Momentum and charge sign of selected neutrino- 
induced muons. The calculated background is shown as well 
as the Monte Carlo expectation in the absence of oscillations. 



Pfit (GeV) 


Data 


Bkgd 


MC 


1 - 10 (L) 




21 


2.2 


37.5 


10-100 (H) 


20 


0.2 


17.5 


1 - 10 (L) 




16 


1.3 


19.3 


10-100 (H) 


13 


0.2 


8.6 


unknown (U) 


U 


70 | 0.7 | 76.5 



A. Systematic Uncertainties 

There are several sources of systematic uncertainty in 
this analysis. These are due to both the event recon- 
struction and the physics modeling. These uncertainties 
are summarized in Table HV1 

The sources of reconstruction systematic uncertain- 
ties are those associated with the data selection cuts 
(1) X%t/ nd f> ( 2 ) xl/p/ndf and (3) x\ nc /ndf, where the 
numbers refer to Table IIV1 The systematic uncertain- 
ties for these cuts were all computed in a similar man- 
ner. For example, we computed the systematic uncer- 
tainty on Xi/p/ndf by establishing that the cut value of 
xl/p/ndf = 3 selects 98.9% of events seen in Fig. [U We 
divided the total Monte Carlo data set into 12 subsam- 
ples of 20,000 events each. For each subsample we found 
the value of xl/p/ndf that accepts 98.9% of the events. 
The variance in these values of x\/pl n df divided by the 
nominal value of the Xi/p/ n df cut is the la uncertainty 
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TABLE IV: Sources of systematic uncertainty in low to high 
momentum event ratio 1Z. 



Source 


a 


Aii k 


Reconstruction systematics: 






(1) Xm/ndf < 1.5 


0.01 


< 5 x 10~ 4 


(2) xt/p/ndf < 3.0 


0.01 


< 1 x 10" 4 


(3) Xime/ndf < 10 


0.27 


0.02 


Model systematics: 






(4) Normalization 


0.15 


< 1 x 10 -4 


(5) Spectral Index 


0.03 


0.08 


(6) Cross Section (E < 30 GeV) 


0.07 


0.03 


(7) Cross Section (E > 30 GeV) 


0.02 


0.01 






0.09 



given in Table IIVI The la values for the uncertainties 
in Xfit/ nc tf an( A X\ine/ n< V were computed in a similar 
manner. 

The first source of systematic uncertainty in the 
physics modeling is the (4) overall normalization of the 
calculated neutrino flux. Uncertainties in the primary 
cosmic ray fluxes and the hadronic production are the 
main contributors to the overall uncertainty in the nor- 
malization. Combining these sources give s an uncer- 
tainty in the flux normalization of 15% [35j. There are 
three sources of systematic uncertainty in the ratio of 
the number of muons induced by low energy neutrinos to 
those induced by high energy neutrinos. One contribu- 
tion comes from (5) the uncertainty in the spectral index 
of the atmospheric neutrino energy spectrum. The neu- 
trino flux is proportional to E u with the value of 
the spectral index, 7 = 1.7 ± 0.05 [HI, a 3% uncertainty 
in the spectral index. In addition, the uncertainties in 
the neutrino and anti-neutrino cross sections contribute 
another 7% [3(| to the rate of muons coming from neu- 
trinos with energies < 30 GeV (6) and 2% [10] for muons 
coming from neutrinos with energies > 30 GeV (7). 

B. Low to High and Unknown Momentum Event 
Ratio 

One way to look for evidence of neutrino oscillations 
in the neutrino-induced muons is to take the ratio of the 
number of low momentum muons, which are more likely 
to show an oscillation signal, to the sum of the number of 
high momentum and unknown momentum muons, which 
are less likely, and compare this ratio with its Monte 
Carlo expectation including backgrounds. In the data, 
this ratio of low to the sum of high and unknown mo- 
mentum muons is given by 

Rt a /H + u = £(^v + jvv E ( N »- + ^v), (2) 

L H+U 

where N„- is the number of [i~ observed in a bin and 
iV^+ is the number of /Lt + observed in a bin. The sum over 
L includes events in the range 1 < pfu < 10 GeV/c, and 



the sum over (H + U) includes the remaining high mo- 
mentum and unknown momentum events. In the Monte 
Carlo simulation, a similar ratio, R^/h+w lS defined. In 
the absence of oscillations, the ratio of these two quan- 
tities, 7Z, will be consistent with unity; if an oscillation 
signal is present, 1Z will be less than unity. 

We computed the systematic uncertainty in this ratio 
with our Monte Carlo simulation by varying the value 
of each of the parameters in Table UVl by ±lcr from the 
values used in the analysis. The change in the Monte 
Carlo ratio resulting from the variation in each parame- 
ter, AlZk, was assumed to be the uncertainty in the ratio. 
The total uncertainty in A1Z was found by adding the un- 
certainties from the individual parameters in quadrature. 

As the reconstruction uncertainties are based on how 
different the cut efficiencies are for data and Monte Carlo 
simulation, we determined how the value of those se- 
lection criteria affected the ratio by varying the value 
of the cuts. Varying j& it /ndf between 1.485 and 
1.515 gives ATZi < 5 x 10~ 4 . The value A1Z 2 < 
1 x 10~ 4 is as expected from Fig. |4j Changing the 
value of Xii nc / n df between 7.3 and 12.7 shows that 
AIZ3 = ±0.02. The flux normalization simply scales 
the number of low momentum, high momentum, and 
unknown momentum muons by a constant so this un- 
certainty cancels in the ratio. Varying the spectral index 
by ±3% alters the relative numbers of low momentum to 
high momentum and unknown momentum muons which 
leads to AIZ5 = ±0.08. Varying the cross section for neu- 
trinos with E < 30 GeV by ±7% gives ATZ 6 = ±0.03. A 
variation of ±2% in the cross section for neutrino inter- 
actions with E > 30 GeV gives A1Z 7 = ±0.01. Adding 
these uncertainties in quadrature gives the total system- 
atic uncertainty, A1Z — 0.09. These results are given in 
Table 13 

From Table [HI] and Table [IV] we find 

j^data 

71 = mc +U = 0.65±g;^|(stat) ± 0.09(syst). (3) 

R-L/H+U 

The upper and lower limits on the data event rate ratio 
are estimated accounting for the Poisson fluctuations in 
the numerator and demoninator [37| . The ranges quoted 
are calculated to give coverage at 68% C. L. Adding the 
upper statistical and systematic uncertainties in quadra- 
ture, the upper uncertainty is +0.17 which results in a 
value for 1Z that differs from the no oscillation expec- 
tation of unity by 2. Oct. This result is consistent with 
neutrino oscillations. 



C. Oscillation Fit 

In the following section we test the significance of the 
neutrino disappearance suggested by the value of 1Z in 
Eq. Q by fitting for the oscillation parameters sin 2 28 2 3 
and Ami. 
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X 2 /ndf = 5.9/7 

MINOS Beam Best Fit 
No Oscillations 
Background 



FIG. 14: The intensity of neutrino-induced muons as a func- 
tion of cos#. The data are shown by the points, the best fit 
is shown by the solid line, and the null oscillation hypoth- 
esis is shown by the dotted line. The prediction using the 
MINOS result with the NuMI neutrino beam is shown by 
the dashed line. The top left panel shows the events with 
1 < Pfit < 10 GeV/c (L), the top right shows the events with 
10 < Pfu < 100 GeV/c (H) and the bottom left shows events 
with unknown momentum (U). 



1. Method 

We used the data and Monte Carlo samples from 
Table IIIII for the oscillation fit. We first divided the 
momentum-separated muons (L,H,U) into three cos# 
bins: -1.0 < cos6> < -0.6, -0.6 < cos6> < -0.2 and 
—0.2 < cos# < 0.05. To find the oscillation parame- 
ters (sin 2 2#23, Am 2 3 ), we compared the data and Monte 
Carlo simulation with the x 2 statistic, 



i=i 



MO 



TV, D + Nj In 



N MC 



fc=i 



(4) 



The first sum is the \ 2 for Poisson distributed data [38| . 
The number of data events seen in bin i is NP and the 
number of events predicted in bin i by the Monte Carlo 
simulation is N^ . The parameters in the second 
sum account for the systematic uncertainties described 
in HIV Al The are included in the fit as penalty terms 
and allow the fit to adjust the predicted number of muons 
in each bin i based on the systematic uncertainties. The 
ek are assumed to be Gaussian distributed with a mean of 
so that excursions of the fit values from the expected 
values contribute an amount e 2 /<7 2 to the x 2 statistic, 
where the cr^ are given in Table HVl 

The iVj MC depend on the effects of neutrino oscilla- 
tions, the systematic uncertainties in the reconstruction 
and input physics model, and the background present for 
bin i. We define the expected number of muons in bin 
* as Nf c = iVf c (sin 2 26> 23 , Am 2 3 , e) + N® G . The back- 



10- 



> 

eg cm 

E 

< 

1( 



10" ! 



k 
; 

^ v 

^ Best Fit 
68%v-induced 

— 90% v-induced 

— - 90% contained 



0.2 



0.4 



0.6 0.8 



FIG. 15: The the 68% (dotted line) and 90% (solid line) con- 
fidence intervals for the oscillation parameter fit. The best 
fit point is indicated by the star. Also shown is the 90% 
confidence interval for the MINOS contained vertex analysis 
(dashed line) based on the first 418 days of data taking with 
the far detector. 



ground events due to neutrino interactions inside the de- 
tector, N® G , have been oscillated with the best fit Super- 
Kamiokande parameters and e = (ei, . . . , £7) correspond 
to the uncertainties listed in Table llVl We calculated the 



at each point to find the 



values of x 2 a t each point on a grid in (sin 2023, Am 2 3 ) 
space. The fit varied the e 
minimum \ 2 [3^ | in Eq. j4|) 



2. Results of Oscillation Fit 



The 
^ 2 28, 



best fit in the physical region is found at 
sin" 2V 2 3 = 1 and Am| 3 = 0.93 x 10~ 3 eV 2 with x 2 / n df = 
5.9/7. The values of the were all found to be less than 
lcrfc. The results of the fit are shown in Fig.[T21 where the 
best fit to the data in the physical region (solid line), the 
Monte Carlo prediction for the null oscillation hypothesis 
(dotted line), and the prediction derived from the MINOS 
result with the NuMI neutrino beam (dashed line) [l9| 
are superposed onto the observed intensity of neutrino- 
induced muons. For the results from the null oscillation 
hypothesis and the MINOS beam analysis we varied the 
ej; to minimize Eq. ^ at the fixed values of sin 2 2#23 and 
Ami,. 

Fig. [TSI shows the 68% and 90% confidence intervals for 
the oscillation parameter fit. The points within the 68% 
contour have a difference in \ 2 with the best fit point, 
Ax 2 , of less than 2.3. The points within the 90% contour 
have Ax 2 < 4.61. The null oscillation hypothesis has 
Ax 2 = 5.5 and is excluded at the 94% confidence level. 
Also shown in Fig. [15] is the 90% confidence interval for 
the MINOS contained vertex analysis based on the 
first 418 days of data taking with the far detector. The 
results of these two analyses are consistent. 
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V. TESTS OF CPT CONSERVATION 

To probe CPT conservation with MINOS data, we use 
two tests. After discussing the systematic uncertainties, 
we first compare the strength of the oscillation signals 
for the neutrino-induced and [i + samples individu- 
ally. Second, we divide the charge-separated muons into 
low and high momentum samples and further test for 
differing rates of neutrino disappearance. We show that 
this second test is the one that minimizes systematic un- 
certainties. 



A. Systematic Uncertainties 

Here we discuss the additional sources of systematic 
uncertainties associated with the CPT analysis. One 
additional source of systematic uncertainty is due to 
misidentified charge sign or the charge purity of the 
data sample. We determined the charge purity system- 
atic by comparing the underground cosmic ray muon 
charge ratio for events passing the Xu nc / n df cut with the 
charge ratio found by MINOS, = 1.371 [H. 

For this comparison we divided the cosmic ray muons 
into two data samples, one with muons having momenta 
< 30 GeV/c and the other with momenta in the range 
30 - 100 GeV/c. For cosmic ray muons in the lower 
momentum sample, we calculated the difference in the 
charge ratio to be 0.047, a difference that can be at- 
tributed to an impurity of 5%. For cosmic ray muons in 
the high momentum sample, we calculated the difference 
in the charge ratio to be 0.12, which can be attributed 
to an impurity of 14%. By taking a weighted average 
of these impurities, where the weights are given by the 
number of neutrino-induced muons in each momentum 
range, we calculated the impurity in the neutrino-induced 
muon data to be 6%. Using similar methods we found 
the weighted average impurity for the ncutino-induced 
Monte Carlo events to be 1%. These weighted impurities 
are given in Table [V] (4) . 

There are also additional sources of systematic uncer- 
tainty due to the physics model. These arise from uncer- 
tainties in the ratios of the cross section, av „ , to the 
cross section, <t„ , and we estimate the uncertainties 
to be 9% [36| for energies < 30 GeV (9) and 2% [U for 
energies > 30 GeV (10). The uncertainty in the ratio of 
the number of to (11) is 1% [35j . 

B. Charge-Separated Event Ratio 

Our first test of CPT conservation comes from the com- 
parison of the ratio of the total number of \x~ events to 
/i + events in the data and the Monte Carlo simulation 
including backgrounds. In the data, this ratio of charge 
separated muons is given by 

^-/V = E(^-)/E( jV ^)- ( 5 ) 



TABLE V: Sources of systematic uncertainty considered and 
their effects on the ratio IZcpt- 



Source 


a 


(A1Z CPT ) k 


Reconstruction systematics: 




(1) Xm/ndf < 1.5 


0.01 


< 5 x 10" 4 


(2) xt/g/ndf < 3.0 


0.01 


< 1 x 10~ 4 


(3) xLJndf < 10 


0.27 


< 5 x 10 _a 


(4) Charge Purity (Data,MC) 


(0.06,0.01) 


+0.07 


Model systematics: 




(5) Normalization 


0.15 


< 1 x 10" 4 


(fit Slnppfrfll TtiHpy 


0.03 


< 1 x 10~ 4 


(7) Cross Section (E < 30 GeV) 


0.07 


< 5 x W- 5 


(8) Cross Section (E > 30 GeV) 


0.02 


< 5 x 10-* 


(9)<W^„ (E<30 GeV) 


0.09 


0.04 


(10) a v Ja Vu (£>30 GeV) 


0.02 


< 5 x 1CT 3 


(11) Vp/Vn 


0.01 


0.01 


ATZcpt = \fj2( A ^-CPT)i 




-0.04, +0.08 



The Monte Carlo ratio is defined in a similar way. If 
the f M oscillate with the same parameters as the V ^ then 
the ratio of the data and Monte Carlo ratios, IZcpt, 
will be consistent with unity. The systematic uncertain- 
ties in licPT are shown in Table fVl This table shows 
how changes in the default parameters corresponding to 
the sources of systematic uncertainty changes the value 
of IZcpt- The main contributions to the uncertainty 
in IZcpt are the purity, cross section ratio for neutri- 
nos with energies < 30 GeV, and the ratio of v^/V^. 
Since the ratio in Eq. (J5]) is greater than one, more nega- 
tive than positive muons will be misidentified and the 
measured charge ratio decreases towards unity. Con- 
sequently, charge misidentification leads to a one-sided 
(positive) error on the measured charge ratio. The un- 
certainty due to the impurity is (A1Zcpt)a = 0.07. The 
cross section ratio for neutrinos with energies < 30 GeV 
changes the ratio by (A1Zcpt)v — ±0.04. The uncer- 
tainty in the ratio of the number of to corresponds 
to (ATZcpt)h = ±0.01. We find from Table [nT] and 
Table E] 

j^data 

Kcpt = = 0.72:^(stat)I°^(syst), (6) 

where the statistical uncertainties were calculated using 
the method for Eq. ([3]) . This value of IZc pt is consistent 
with CPT conservation. 



C. Charge-Separated, Low-to-High Momentun 
Event Ratio 

Our second test of CPT conservation is based 
on the charge separated ratios of low momentum to 
high momentum events compared to the Monte Carlo 
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expectation including backgrounds. First define ratios 
of data to Monte Carlo expectation for the and 
individually, 



jjdata 
L/H 

n L/H 



0.49 



-0.20 
-0.14 



(stat) ±0.07(syst), 



(J) 



and 



r>data 
L/H 

L/H 

= 0.55±g'jg(stat) ± 0.07(syst). 



(8) 



The systematic uncertainties in these ratios are listed 
in Table IVI| and the statistical uncertainties were again 
calculated using the same method as for Eq. ([3]). We 
define the ratio of ratios IZcpt, 

Kcpt = 11-/11+ = 0.89±g;^(stat) ± 0.03(syst). (9) 

We found the upper and lower limits for the statistical 
uncertainty in Eq. ^ by fitting the observed L and H 
event rates for ^~ and [i + to a model which used IZcpt 
and three of the four event rates as inputs. The data and 
model were compared using a x 2 function appropriate 
for Poisson distributed data (38|. The 68% C. L. interval 
quoted includes all values of IZcpt which produce a x 2 
within 1 of the best fit value. Again within the statis- 
tical and systematic uncertainties, the value of IZc pt is 
consistent with unity and CPT conservation. 

The value of this test is clearly demonstrated in Ta- 
ble IVII where it can be seen that several sources of sys- 
tematic uncertainty cancel in the calculation of IZcpt- 
These cancellations occur because the fi~ and /i + ratios 
are affected similarly for these systematic uncertainties 
and so divide out in the ratio. The uncertainty in the 
ratio of the number of to the number of does not 
contribute to the uncertainty in IZcpt- As seen in Ta- 
ble |Vll the only sources of uncertainty which do not can- 
cel are those in the relative magnitude of the and 
cross sections at low (9) and high (10) energies. 

This approach should prove valuable in high statistics 
accelerator-based tests of CPT conservation because it 
minimizes systematic uncertainties so effectively. 



VI. CONCLUSIONS 



The MINOS far detector has taken 854.24 days of data 
in its search for the atmospheric neutrino- induced muons. 
We found a total of 140 neutrino-induced muons in the 
data set. We have performed four analyses. The ra- 
tio of the number of low momentum muons to the sum 
of high momentum and unknown momentum muons has 
been compared to the same ratio for Monte Carlo events. 
The resulting ratio, 1Z, suggests neutrino disappearance 
at the 2.0cr level. The best oscillation fit to the data gives 
sin 2 2023 = 1.0 and Am| = 0.93 x 10~ 3 eV 2 . This fit 
excludes the null oscillation hypothesis at the 94% con- 
fidence level. We used two methods to look for evidence 
of different rates of neutrino disappearance as a test of 
CPT conservation. Both tests are consistent with CPT 
conservation. The cancellation of many systematic un- 
certainties in the computation of IZc pt suggests that this 
test would provide a precision test of CPT conservation 
with a suitably large dataset. 
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